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Abstract 

It is shown that in 2D system of the fermions with simplest indirect boson-induced 
attraction (through the Einstein phonon exchange as an example) along with the nor- 
mal and superconducting phases there arises a new (called "abnormal normal" or pseu- 
dogap) one where the absolute value of the order parameter is finite but its phase is 
a random quantity. It is important that this new phase really exists at low carriers 
density only, i.e. it shrinks when doping increases. The relevance of the results ob- 
tained with dependence of pseudogap on doping in high-temperature superconductors 
is speculated. 
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1 Introduction 



An adequate description of the physical properties of high-temperature superconductors 
(HTSCs) is one of the most important problems of the modern solid state physics (see, for 
example, review [![]). And among the most debatable questions on HTSCs is the question 
about the so called pseudogap (or spin gap) which is experimentally observed in normal 
state samples with lowered carrier density rif. The matter is that HTSCs represent systems 
with relatively easy changeable value of rif what is interesting and important both from 
experimental and from theoretical points of view. In particular, the underdoped samples 
reveal a lot of strange anomalies concerning their spectral, magnetic, thermodynamic and 
other observables (among them there is the pseudogap which has been seen in ARPES 
experiments |2|, [| |J] and now widely discussed || [£], 0). Moreover, even the definition 
itself "underdoped" ("overdoped") HTSC sample is related to the presence (absence) of 
such anomalies. 

On the other hand, the possibility of changing rif value is crucial for the theory and 
puts the general problem of the crossover from composite boson superfluidity (underdoped 
regime, or small rif) to Cooper pairing when rif increases (overdoped regime). This crossover 
was already studied for the systems of different dimensionalities: 3D || and quasi-2D ||. As 



for 2D case, it has been considered for temperature T = only || |T0| what is conditioned by 
the Hohenberg-Mermin- Wagner theorem which forbids any long-range order in such systems 
at finite T because of long-wave fluctuations of the charged order parameter. 

Formation of inhomogeneous condensate, or the Berezinskii-Kosterlitz-Thouless (BKT) 
phase, is possible, but its consistent study is connected with some difficulties. Nevertheless, 



this case was already explored in 2+1 relativistic field model |TT| where, as it is known, 
concentration effects are not justifiable. They are known for some nonrelativistic models 
(see, for example, |12|]), but the BKT equation in this paper was obtained without taking 
into account the existence of neutral order parameter p 0). Its introduction and use, as 



it was shown in [I3"|, is very important and leads to formation of a new phase with p ^ 
which separates normal phase and another one, which is also normal, since superconducting 
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properties are absent there as well. 

Below an attempt is made to study the above mentioned crossover and the new phase 
formation possibility. In contrast to |T3, where this question was studied for a 2D four- 



fermion (4F) model, we shall consider a more realistic Frohlich model. It will be shown that 
the unknown phase in this case not just exists, but appears mainly at low carrier density 
what is essentially different from the 4F-case. Such a result promises that the models with 
an indirect attractive femion-fermion interaction may be suitable to account for the unusual 
features in the normal phase behaviour of HTSCs. 

2 Model and main equations 

Let us choose the Frohlich model Hamiltonian density in the standard form: 

H(x) = -tpl(x) f ^- - [A tj) a {x) + gip(x)ipl(x)ip a (x) + H ph (ip(x)), (x = r,t), (1) 

where ip a (x) is a fermion field; m and a =f, 1 are an effective mass and a spin of the fermions, 
H is their chemical potential which fixes <£>(x) is a phonon field operator and g is a fermion- 
phonon coupling constant. We use Pauli matrices J, r x , r y , t z and put % = ks = 1. 

In ([!]) H p h is the Hamiltonian of free phonons, which results in the simplest expression 
for the phonon propagator (in the Matsubara temperature formalism) [[14]] 

to 2 

D(iV n ) = - — 2 , n n = 2nnT (2) 

with Uq as Einstein phonon frequency and n as an integer. In principle, @ corresponds to 
a propagator of any massive bosonic excitations by which fermions can exchange^. 

With the purpose to calculate the phase diagram it is necessary to find the thermody- 
namic potential of the system. In the case of the 4F model it can be obtained by making use of 
1 In spite of its simplicity such a model is rather close to HTSCs with their developed quadrupolar mode 
spectra - optical phonons and dcZ-excitons Q . Exchange of these excitations leads to short-range interaction 
between fermions. It, however, does not include spin excitations (almost massless in HTSCs) which can 
provide the long-range inter-carrier interaction what, of course, needs a separate study. 
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the well-developed Hubbard-Stratonovich method in which the statistical sum Z can be rep- 
resented by a path integral over the fermionic ip^{x) and the complex Hubbard-Stratonovich 
4>{x) = V < Vwl > fields (V is the 4F-interaction constant). 

As it was shown earlier [p~|, [L3| , in 2D case it is convenient to pass from <p and (f>* fields to 
new variables, namely: the absolute value p and the phase 9, where <j>{x) = p(x) exp[— 2i9(x)], 
and to perform simultaneously the spinor transformation 

^a(x) = Xa(x)exp[i9(x)], (3) 

where x<r{ x ) corresponds to a neutral fermi-particle field. Such a substitution allows to 
represent Z in the following form 



pVpVO exp 



p / tt(p(x),6(x))dx 



(P = (4) 



where Q is the effective quantum thermodynamic potential in terms of p(x) and 8(x) vari- 
ables. This representation proves to be very useful and leads to (at p(x) = const and 
expansion in V6 l (x) up to (V6(x)) 2 ) the effective Hamiltonian which is similar to that of the 
XY-model. As a result, the equation for Tbkt i n the latter model can be directly used for 
the case under consideration. The solution of whole the set of the self-consistent equations 
for p, p and T BKT was given in [ |3~3f1 . 

This approach is however inapplicable in the more complex case of indirect interac- 
tion model for which the local Hubbard-Stratonovich fields can not be introduced. At the 
same time, by making use of the Cornwall- Jackiw-Tomboulis formalism [Tj| it is possible 
to obtain classical thermodynamic potential which depends on bilocal < ipi(x)ip^(y) > and 
< il)^{x)ip\{y) > fields (or, more exactly, on the full fermion Green function (see below)) 
though, strictly speaking, this potential does not satisfy Eq. (J|) exactly. The corresponding 
effective action (equivalent to Q in (f|)) in two- loop approximation can be calculated and (in 
the Nambu representation |17|) takes the form 



(3Q[g] = -Tr (Lng-'Qo + QQ^ - lj + »-^Qt z Dt z Q, (5) 
where Tr includes integration over 2D space r and imaginary time < r < ft as well as the 



standard trace operation; 



£-! = -W t + t z f^ + /i) + r + 0* + r_0; (6) 



2m 



g-i = g-\ P = o) 



are the full and free fermion Green functions, and D was defined in (|2|). We also use the 
normalization condition Q = Qo(fi = 0) under the symbol Ln in (H). 

The stationary condition 5Q(Q)/5Q = results in the well-known Schwinger-Dyson equa- 
tion for the inverse full fermion Green function 

Q- 1 = Go 1 ~ 9 2 TzGr z D. (7) 

Substituting (171) into (0) one can obtain a simpler expression for where the phonon 

Green function (|2]) is already omitted 



(3Q(g) = -Tr 



Lner 1 ^ + - [gg- 1 - 1 



To find the potential Q as a function of phase gradient V6(x) and absolute value p it is 
necessary to make the transformation @. Then one obtains for (^): 

-1 /V 2 \ 

g = -Id T + r-A — + (J, I + r x p - 

So" 1 = G^-EW 

Since the low energy dynamics is mainly determined by the phase fluctuations and cor- 
responds to the region where spatially homogeneous order parameter p ^ it is sufficient 
to be restricted the expansion in terms of Q in V8(x) only; for instance, the Green function 
([]) g = G + G J2^Lo(^^) n ■ Then the desirable effective potential (^|) can be divided it two 
parts: Q = Q pot (p) + Q kin (p, V#), where in (V6>) 2 approximation 



pn kin (p,V6) = Tr 
P 



G£ + ^GSGS - G S - ^GoSG S+ 



t x G(GT, + GSGE) 



(10) 



2 In fact it coincides with the Eliashberg equation |l7|] when the renormalization factor is Z(iu> n ) = 1. 



Assuming then that p{iuj n ) = const one obtains from ( |T0D after somewhat tedious but 
otherwise straightforward calculation 

T rP 
2 

where 



tt kin = \\l dr \ d 2 rJ(p,T,p(p,T))(V9f 



J(p,T,p(p,T)) = — U 



T 
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2 + p 2 + p + 2T\n 

2 Q 



P 



cxp 



dx- 



x + (n/2T) 



M/2T cos h 2 J x 2 + p 2 /AT 2 



(12) 



AT 2 d(p 2 /AT 2 ) 

Note that in comparison with the 4F model the last expression contains one more term 
with the derivative. 

The equation for the temperature Tbkt of the BKT transition can be written after direct 
comparison of Qkin (0) with the Hamiltonian of the XY-model which has the same form 
HH; hence 

71 

-J{p, Tbkt, p(p, Tbkt)) = T B kt- (13) 

To complete a set of selfconsitent equations which allow us to trace the dependence of Tbkt 
on rif the equations for p and p have also to be given. 

In particular the equation for p(iou n ) is nothing else but (0) with V0 = 0, i.e. the 
Green function G of the neutral fermions replaces Q, so that ([7]) in frequency-momentum 
represantation takes the form 



00 . 

p{iu n ) = g 2 T J 



d 2 k 



p(zw r , 



(2vr) 2 Lo 2 m + e 2 (k) + p 2 {iu m ) (uj m - u n ) 2 + u% ' 



(14) 



where u n = (2n + 1)ttT is the Matsubara fermionic frequency JT] 

Analytical study of this equation as well as obtaining both ( |T2"D and the number equation 
is possible only if one supposes that p{iu) n ) = const 

Finally, the number equation which follows from the condition v~ 1 dQ[Q]/dp = —nf 
(where v ia a volume of the system) and is crucial for crossover description has to be added 
to ( |i~3"D and fll4]) for self-consistency; so one comes to 

VJ^+p 1 ' 



p 2 + p 2 + p + IT In 



1 + exp 



T 



2e r 



(15) 



3 For the case T = this assumption was checked in |lj 



where ep = irnf/mis the Fermi energy of 2D fermions with the simplest quadratic dispersion 
law. Thus, in the case under consideration all unknown quantities (p, p and Tbkt) are 
functions of rif. 



3 Analysis of the solutions 

It is natural to guess that in contrast to the accepted XY-model in the superconducting 
model there exists one more critical temperature, also dependent on rif, where the complete 
order parameter vanishes. This critical temperature, T p , can be found from flU]) and ( |15|) 
by putting p = (what in accordance with these equations derivation corresponds to the 
mean-field approximation). As a result, with temperature decreasing 2D metal passes from 
normal phase (T > T p ) to another one where averaged homogeneous (charged) order 
parameter < <p(x) >= 0, or, what is the same, superconductivity is absent, but chargeless 
order parameter p ^ 0. It is very important that the pseudogap is formed just in the 
temperature region T BKT < T < T p , because, as it follows from above formulas (see, for 
instance, (0), ([Top), p, which depends on rif and T, enters all spectral characteristics of 
2D metal in the same way as the superconducting gap A(T) enters into corresponding 
expressions for ordinary superconductors. It explains why this, new phase can be called 
the "abnormal normal" phase or, better, pseudogap one. The density of states near ep 
in the pseudogap is definitely less than in the normal phase, but is not equal zero as in 
superconducting one. 

The phase diagram of the system is presented in fig.l. Both curves for T p and Tbkt as 
functions of rif can be obtained by numerical calculation. As for their asymptotic behaviour 
can be established analitically. 

Indeed, it is not difficult to be convinced that the asymptotics for T p (rif) and T BKT {rif) 
have the forms: 

a) when ratio tpjuj^ <C 1 (very low fermion density, or local pair case) the first of them sat- 
isfies the equation T p \n(T p /e F ) = cj exp(— 2/A) which immediately results in dT p (rif) /drif\ nf ^ 
oo. At the same time the temperature Tbkt here has has another carrier density depen- 
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dence: Tbkt = £f/2 what in its turn means that it is equal to the number of composite 
bosons. It means that in this density region T p /Tbkt S> 1- 

b) in the case ep/ujQ ^> 1 (very large fermion density, or Cooper pair case) one comes 
to standard BCS value: T p = {Z^ujq/k) exp(— 1/A) = Tqq S = {2^ /n)A. B cs (^bcs is the 
one-particle BCS gap at T = 0). In other words, the temperature T p becomes equal to its 
BCS value 0. The Tbkt asymptotics in this case is not so evident and requires more detailed 
consideration. 

Firstly, it is naturally to suppose that for large n/ values Tbkt T p . Then it is necessary 



to check the dependence of p on T as T — > T p . For that the equation flT4j ) can be transformed 
to: 

1 f oo ( tanh ^x 2 + p 2 /4T 2 tanh ^x 2 + p 2 /4T 2 - tanh(w /2T) 
dx 1 



A 



V Jx 2 + p 2 /4T 2 2(Jx 2 + p 2 /AT 2 -u /2T) 



tanh Jx 2 + p 2 /AT 2 + tanh(cu /2T) 



(16) 

2(tJx 2 + p 2 /4T 2 + uj /2T) J 
(where it was used that in this concentration region the ratio p/2T p ~ ep/2T p ^> 1). 

Because usually ujq/2T p ^> 1, only very small x give the main contribution to the integral 
(|16|) (it is seen from the limit p/2T p — > when e^/cuo ~^ oo). Therefore the latter expression 
takes the approximate form: 



1 = « tanh^ + p 2 ^ 2 !_ (i7) 

A Jo I Jx 2 + p 2 /4T 2 x + u /2T> 



On the other hand, the condition p = in (|T7D leads to the equation 

1 r°° , ( tanh x 1 \ 

dx 



A Jo \ x x + ujq/ 2T P t 
for T p . From (|T7|) and ( |18D it directly follows that 



o 



/ tanhx tanh^ + p 2 ^ 2 ^ _ i T, 

I I I — 111 • 

V x Jx 2 + p 2 /4T 2 j T 



4 Being equal (in mean field approximation only) these temperatures (T p and Tg^g) are in fact different: if 
Trcs immediately falls down to zero as fluctuations are taken into account, T p does not and is renormalized 
only. 
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Then using the approximation 



tanh Jx 2 + p 2 /AT 2 1 - 3" 1 [x 2 + p 2 /AT 2 } , x < 1 



x 2 



p 2/4T 2 [ x -i - p 2 /8T 2 x 3 , x > 1 

one easily comes to expression needed: 



p~2.62T p JT p /T-l. (19) 



Remind that the generally accepted 3D result is A BC s(T) = 3.06Tj^y Tj^/T — 1 [14] and 
this small difference can be explained by the above approximation what, however, is suitable 
for the following below qualitative discussion (see next Section). 

The last dependence has to be substituted in equation ([T3|). And because of p/2T BKT ~ 
ep/2T BKT ^> 1, and p(T BKT )/2T BKT <C 1 when T BKT — > T p this equation can be written as 



AT, 



BKT 



p 2 d 



±TIkt d(p 2 /AT 2 BKT ) 



dx 



o V cosh 2 x 



1. (20) 



cosh 2 + p 2 / 'AT 2 KTi 
At last, using expansion in p/2T BKT in integral fl20D , the latter can be transformed to 



ae B / p \ 4 



8T BKT \^.i. L;KT 

where the numerical constant 



^— =1, (21) 



f°° , tanh 2 x — x 1 tanh x + 1 
a = I dx — ~ 1.98. 

2x z cosh x 







Combining now (O) and (^1|) one comes to the final relation between T p and T B kt for large 
carrier density: 



T BKT ~T p (l-1.17^/T p /e F ), (22) 
i.e. Tbkt as a function of rif approaches T p (or T^fg). 



4 Discussion 



Thus, in the considered model of 2D metal with indirect interection of carriers, along with the 
normal phase (p = 0), there exist pseudogap and superconducting (here - BKT) phases which 
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in fact correspond to the cases of absent, uncorrelated and correlated pairs, respectively. De- 
spite the simplicity of the model and assumptions made above, the strong dependence of the 
critical temperatures T BKT and T p on carrier density is evidently shown. Pseudogap phase 
really exists when e F < lOOTg^ (or (T p — T BKT )/T p > 0.1 (see (|22|))). This conclusion 
agrees qualitatively with the existence of unusual spectral peculiarities in some underdoped 
HTSCs (for example, YBCO || |[ [§]). There are many speculations about the nature of 
this pseudogap phase ( in particular, the spin-gap conception (see ^(|) is also very popular). 
Our results give hopes that "the low-density phenomena" in HTSCs can be (at least, prin- 
cipally) explained in the framework of a rather standard boson-exchange approach though 
the simplest phonon model considered can hardly be available for description of real HTSC 
compounds. 

It must be also noted that unlike this model the 4F one predicts the pseudogap phase 



existence at any fermion densities (T^|. It means that in the case 2D boson-exchange model 
with, as shown above, "saturated" critical temperatures as functions of doping, the defini- 
tions of "underdoped" and "overdoped" samples aquire their physical sense. Nevertheless, 
a point where both (T p (rif) and 7rr-t(^</)) curves cross each other is absent. The matter is 
that the second temperature appears at finite value of p only. In this connection it is worth 
to mention that real HTSCs are not pure 2D, but quasi-2D systems with very weak (mainly, 
Josephson) tunneling between neighbouring conducting copper-oxide layers. In such a case 
there can exist one more critical temperature || 

T „ Tbkt (03) 

Ic -He F \e b \/4qy {Z6) 

which (at small nf) defines the formation of usual homogeneous condensate where not only 
p but also the phase 9 does not depend on coordinates. In (p3|) t\\ is the inter-plane hopping 
(tunneling) constant, and Sb is the two-fermion bound state energy |], ||. When T c < Tbkt 
all three temperatures have physical meaning, but when T c > Tbkt (or, as it follows from 
(fZ3|), t\\ > 2Abcs) the BKT phase cannot be formed and pseudogap phase at T = T c 
transfers directly to the ordinary superconducting phase. 

It cannot be excluded that just such a scenario takes place in HTSC copper oxides. Be- 
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sides, superconducting phase which appears at T = T c does not need a finite p, so the curves 
T p (nf) and T c (nj) can in principle cross each other in some point rif = n°f Pt which, thus, sep- 
arates two different doping regions. One more result following from such a consideration is 
that the ratio 2A/T C proves to be decreasing function of rif approaching this ratio canonical 
BCS value 3.52 at large rif only what agrees with experimental observation (see f2T| ). 

Although above it was pointed out that the pseudogap phase - BKT phase transition 
takes place when p ^ and fluctuations of p are not so important as those of 8, it would be 
interesting to estimate their contributions separately. The one-particle gapless spectrum of 
the normal phase and especially pseudogap in abnormal normal phase spectrum are also of 
great interest and need further investigation. 
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Figure 1: Phase T -nf diagram of 2D metal with indirect inter-carrier attraction. The solid 
lines correspond to the calculated (for A = 0.5) parts of the functions T p (uf) and 7Wt(?v), 
dashed ones guide on the eye. I, II and III show the regions of the normal, abnormal normal 
and superconducting phases, respectively. 
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